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ABSTRACT 

The robust estimation of the tiny distortions (shears) of galaxy shapes caused by 
weak gravitational lensing in the presence of much larger shape distortions due to the 
point-spread function (PSF) has been widely investigated. One major problem is that 
most galaxy shape measurement methods are subject to bias due to pixel noise in 
the images (“noise bias”). Noise bias is usually characterized using uncorrelated noise 
fields; however, real images typically have low-level noise correlations due to galaxies 
below the detection threshold, and some types of image processing can induce further 
noise correlations. We investigate the effective detection significance and its impact 
on noise bias in the presence of correlated noise for one method of galaxy shape 
estimation. For a fixed noise variance, the biases in galaxy shape estimates can differ 
substantially for uncorrelated versus correlated noise. However, use of an estimate 
of detection significance that accounts for the noise correlations can almost entirely 
remove these differences, leading to consistent values of noise bias as a function of 
detection significance for correlated and uncorrelated noise. We confirm the robustness 
of this finding to properties of the galaxy, the PSF, and the noise field, and quantify 
the impact of anisotropy in the noise correlations. Our results highlight the importance 
of understanding the pixel noise model and its impact on detection signihcances when 
correcting for noise bias on weak lensing. 

Key words: data analysis - techniques: image processing - gravitational lensing: 
weak 


1 INTRODUCTION 

Gravitational lensing, the deflection of light by mass, 
has been used for many different applications due 
to its sensitivity to all gravitating mass, including 
dark matter. Weak gravitational lensing (for a re¬ 
view, see Bartelmann & Schneider 2001; Refregier 2003; 
Schneider 2006; Hoekstra & Jain 2008; Massey et al. 2010a; 
Weinberg et al. 2013), which causes tiny but coherent dis¬ 
tortions in galaxy shapes, can reveal the dark matter halos in 
which galaxies live (e.g., Leauthaud et al. 2012; Tinker et al. 
2013; Velander et al. 2014; Coupon et al. 2015; Han et al. 
2015; Hudson et al. 2015; Zu & Mandelbaum 2015), can 
constrain the amplitude of matter fluctuations and re¬ 
veal the impact of dark energy on structure growth (e.g., 
Hey mans et al. 2013; Jee et al. 2013; Mandelbaum et al. 
2013), and teach us about the theory of gravity on cosmo¬ 
logical scales (e.g., Reyes et al. 2010; Simpson et al. 2013; 
Pullen et al. 2015). 

As lensing surveys grow and their statistical errors de- 
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crease in size, systematic errors become more important. 
Among the main systematic errors in weak lensing measure¬ 
ments are those that arise when trying to accurately infer 
the shears (shape distortions) based on galaxy shape mea¬ 
surements (e.g., Heymans et al. 2006; Massey et al. 2007; 
Bridle et al. 2010; Kitching et al. 2012; Mandelbaum et al. 
2015). An important component of this problem is the sub¬ 
stantial bias that arises due to pixel noise (“noise bias”) 
(Hirata et al. 2004; Kaeprzak et al. 2012; Melchior & Viola 
2012; Refregier et al. 2012), which is typically larger than 
the statistical uncertainties expected from weak lensing 
datasets in the next few years. 

The noise model in astronomical images is typically a 
combination of Poisson noise on the pixel counts plus other 
sources of noise that are Gaussian (e.g., read noise). The 
noise can commonly be modeled as stationary on the scale 
of galaxy images, at least for the faint galaxies that dominate 
weak lensing measurements. In this regime, the variance is 
the same in each pixel, which is true when the Poisson noise 
on the sky level dominates and when the sky level does 
not vary much across each galaxy. In an unprocessed im¬ 
age, the noise is largely uncorrelated between pixels, though 
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the light profiles of galaxies that are just below the detec¬ 
tion threshold do induce some noise correlations. Even more 
significant correlations between the noise in different pixels 
can be induced via processes such as correction for charge 
transfer inefficiency (Massey et al. 2010b) and image resam¬ 
pling (Lauer 1999; Fruchter & Hook 2002; Fruchter 2011; 
Rowe et al. 2011) as part of the image combination process, 
particularly if the images are resampled to smaller pixels 
than those on the original detector. In combined images from 
the HST (Hubble Space Telescope) used for weak leasing, 
this effect can be quite significant (e.g., Koekemoer et al. 
2007), enough so that naive estimates of signal-to-noise ra¬ 
tios are wrong by factors of more than two if the correlations 
are ignored (Casertano et al. 2000). 

In this paper, we therefore consider the issue of noise 
bias in the context of correlated noise. There are several (es¬ 
sentially interchangeable) ways to think about noise bias. In 
general, it can be thought of as a change in shape of the like¬ 
lihood surface due to noise (Refregier et al. 2012), such that 
a maximum likelihood estimator of the per-object shape is 
a biased estimator. Alternatively, for the case of moments- 
based estimators, the bias arises because PSF correction in¬ 
volves the division of two noisy quantities, which has an 
expectation value that is not mathematically equivalent to 
the ratio of their (noiseless) expectation values. The size of 
the bias can be written as a Taylor expansion in the inverse 
detection significance (e.g., Hirata et al. 2004). 

In light of this background, the following questions are 
of particular interest. First, does noise bias have any fun¬ 
damentally different characteristics in the presence of cor¬ 
related noise? Given the mathematical origin of the effect, 
it seems that the noise bias should be described the same 
way for uncorrelated and correlated noise, provided that the 
detection significance is properly quantified. Given this re¬ 
quirement, the second question is whether one can use sim¬ 
ulations with uncorrelated noise to calibrate noise bias even 
for data with correlated noise after accounting for the way 
correlated noise modifies the signal-to-noise ratio (SNR) of 
the object detections? Third, should such corrections use the 
SNR of the galaxy size, flux, or shape? Finally, if correlated 
noise has some directionality, how large is the expected bias 
(Mandelbaum et al. 2012) in shear estimates? We address 
these questions using GalSim (Rowe et al. 2015), an open- 
source^ image simulation package that was designed for tests 
of weak leasing measurement algorithms. It includes rou¬ 
tines to simulate data with uncorrelated or correlated noise 
with a particular correlation function, making it ideal for 
this application. 

In case the answers to these questions depend on the 
galaxy shape measurement method, we use two quite dif¬ 
ferent methods. The first is a moment-based routines (re- 
Gaussianization; Hirata & Seljak 2003) used for the major¬ 
ity of this work, but a subset of our measurements use a 
forward model-fitting method, imSshape (Zuntz et al. 2013, 
2014). A comparison of results with these methods may give 
some insight into the general applicability of our results. 

The outline of this paper is as follows. Sec. 2 defines the 
quantities discussed in the rest of the paper. Sec. 3 includes 
a description of the simulations used for the noise bias in¬ 
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vestigations. Our results are in Sec. 4, with a discussion of 
the implications for a real data analysis in Sec. 5. 


2 TECHNICAL BACKGROUND 

In this section, we define the quantities that are commonly 
measured from galaxy images used for weak leasing, along 
with the uncertainty on those quantities. 


2.1 Galaxy properties 

2.1.1 Galaxy shape 

The majority of the methods of estimating weak leasing 
shear distortions from imaging data involve a two-step pro¬ 
cess (see, e.g., Mandelbaum et al. 2015). The first step is 
to estimate a per-galaxy shape (defined in one of several 
ways), and the second step is to combine all the per-galaxy 
shapes to estimate the cosmological shear or its correlation 
function for the ensemble. The prevalence of the first step 
in methods that have been used for realistic weak leasing 
analysis^ motivates our focus on per-galaxy shape estimates 
as an intermediate step towards the end goal of inferring 
leasing shear. 

Weak gravitational leasing can be described as a linear 
transformation between unlensed and leased coordinates, as 
encoded in the two components of the complex-valued leas¬ 
ing shear 7 = 71 + 172 and the leasing convergence k. The 
shear describes the stretching of galaxy images due to leas¬ 
ing, while the convergence n describes a change in appar¬ 
ent size and flux for leased objects at fixed surface bright¬ 
ness. In practice, the observed quantity is the reduced shear, 
gt = 7i/(l - «)• 

Estimation of shear typically starts with some mea¬ 
sure of the galaxy shape. Despite the fact that galaxies do 
not in general have elliptical isophotes, galaxy light profiles 
are typically modelled as having a well-defined ellipticity, 
£ = {si,S 2 ), with magnitude jej = y/sf s 2 = (1 —5/a)/(l + 
b/a), where b/a is the semi-minor to semi-major axis ratio, 
and orientation angle cj) of the major axis with respect to 
some fixed coordinate system. The two shape components 
can be defined as ei = |£|cos( 20 ) and £2 = |£:| sin ( 20 ), 
or combined as a complex ellipticity s = Si 162 - For a 
randomly-oriented population of source ellipticities, the en¬ 
semble average ellipticity after leasing is an unbiased esti¬ 
mate of the reduced shear: (e) ^ g. 

Another common choice of shape parametrization is 
based on second moments of the galaxy image, 

^ _ J d'^xI{x)W{x)xiXj 
f d2a:7(x)T4^(x) ’ 

where the coordinates xi and X 2 correspond to the x and y 
directions (respectively), I (x) denotes the galaxy image light 
profile, lT(x) is a weighting function (see Schneider 2006), 
and the coordinate origin x = 0 is at the galaxy image center 


^ For alternatives that are under development, see, e.g., 
Bernstein & Armstrong (2014) and Zhang et al. (2015). 
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(the centroid). A second definition of ellipticity, sometimes 
referred to as the distortion, can be written as 


e 


ei + ie2 = 


Qii — Q22 + 2iQi2 
Qii + Q22 


( 2 ) 


For many weight functions W, an image with elliptical 
isophotes of axis ratio b/a has 


_ 1 - 

1 + b'^/a^ ' 


(3) 


For a randomly-oriented population of source distortions, 
the ensemble average e after leasing gives an unbiased es¬ 
timate of approximately twice the shear that depends on 
the population root mean square (RMS) ellipticity, (e) ~ 

2[1 - ((e(^^)^)]s. 

Finally, in the Gaussian approximation 
(Bernstein & Jarvis 2002; Refregier et ah 2012), the statis¬ 
tical uncertainty on e and £ are ae = 2(1 — e^)/z/ ~ 2/z/ 
and as = (1 — £^)/z/ ~ l/u. Here u is a. total detection 
significance. Using simulations with many noise realizations 
to estimate the uncertainty on e allows for the definition of 
a shape-based detection-significance In practice, we use 
the approximate relations above (ue = 2/cre, without the 
factor of 1 — e^) since per-object scatter in e would result 
in very uncertain i/ estimates; however, when using the true 
value of e, the conclusions of the paper are unchanged. 


2.1.2 Galaxy shape estimation methods 

Our work relies on two methods of galaxy shape estimation, 
re-Gaussianization (Hirata & Seljak 2003) and imSshape 
(Zuntz et al. 2013, 2014), that use a distortion e and el¬ 
lipticity s (respectively) to characterize shapes. 

The re-Gaussianization method is a modified version 
of ones that use “adaptive moments” (which are equivalent 
to fitting the light intensity profile to an elliptical Gaus¬ 
sian). The re-Gaussianization method involves determining 
shapes of the PSF-convolved galaxy image based on adap¬ 
tive moments and then correcting the resulting shapes based 
on adaptive moments of the PSF. There are also additional 
steps to correct for non-Gaussianity of both the PSF and the 
galaxy surface brightness profiles (Hirata & Seljak 2003). 

im3shape is a maximum-likelihood fitting method. It 
iteratively fits the selected type of galaxy model to the ob¬ 
served galaxy image, using an oversampled image of the 
PSF. Our adopted settings for all im3shape calculations in 
this work are given in Appendix A. 


2.1.3 Galaxy size 

Many forward-modeling methods of estimating galaxy 
shapes also include an estimate of the intrinsic galaxy size 
(before convolution with the PSF; see, e.g., Zuntz et al. 
2013). Estimates of galaxy intrinsic size are also useful for 
estimates of leasing magnification that take advantage of 
size information (Schmidt et al. 2012; Huff & Graves 2014). 
With im3shape, the natural size estimate to use is the in¬ 
trinsic half-light radius, ri/ 2 - 

While measurements of the observed (PSF-convolved) 
size cr from the elliptical Gaussian-weighted adaptive mo¬ 
ments (Hirata & Seljak 2003) do not say much about the 


real galaxy light profile, they can be useful for quantifying 
the detection significance: 


^size 


a 

Uncertainty on a 


(4) 


2 . 1.4 Flux 


Another type of detection significance is derived from the 
flux from the adaptive moments, and its uncertainty. Note 
that the adaptive moments flux may differ from the total flux 
due to PSF-convolved galaxies not having Gaussian profiles. 
Nonetheless, it is a fast and efficient way to estimate another 
form of detection significance: 


^flux 


flux 

Uncertainty on flux 


(5) 


2.2 SNR estimators 


Estimates of noise bias using simulations are typically done 
as a function of the SNR of the object detection. Sev¬ 
eral weak leasing community challenges (Bridle et al. 2009; 
Mandelbaum et al. 2014) have used an optimal SNR esti¬ 
mator that uses the object light profile itself as a weight 
function, which is only possible if the true noise-free light 
profile is known. This optimal estimator, r'ldeai, is defined as 
follows. The signal S can be defined using a weighted sum 
over the pixels in the image. 


^ Ety(x)/(x) 
EW"(x) 


(6) 


and its variance is 


Var(5) = 


X;vy^(x)Var(7(x)) 

(EW^(x))2 


(7) 


In the limit that the sky background dominates, Var(/(x)) 
is a constant, denoted Var(/(x)) = (the pixel noise vari¬ 
ance). Adopting a matched filter for W, i.e., IU(x) = /(x) 
and putting our assumptions into Eqs. (6) and (7) gives 


^ideal 


(Tp 


( 8 ) 


However, with noisy images, z^ideai is not a measurable 
quantity, forcing us to use realistic estimators of SNR that 
can differ from z^ideai by as much as a factor of two (see, 
e.g., Mandelbaum et al. 2015). As noted in Sec. 2.1, there 
are several possible empirical estimates of detection signif¬ 
icance: z^e, i^size, and z/flux- These can be estimated using 
analytic formulae in some cases, or using measurements of 
multiple noise realizations for each object. One of our goals 
is to compare the behavior of these different estimators of 
detection significance and relate them to z/ideai- 


2.3 Pixel noise models 

The noise fields used in this work are stationary (noise vari¬ 
ance constant across the image) and Gaussian-distributed. 
This approximation is appropriate in the limit that the sky 
background dominates and is large enough that a Poisson 
distribution is nearly Gaussian, which describes most galax¬ 
ies used for weak leasing measurements from the ground and 
a smaller but non-negligible fraction from space. 
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Our uncorrelated noise fields are simply drawn from 
the appropriate Gaussian distribution, with a new random 
number drawn for each pixel. The correlated noise fields 
include the full scale-dependent and direction-dependent 
correlation between noises in adjacent pixels in the im¬ 
ages used for weak leasing science in the COSMOS sur¬ 
vey (Koekemoer et al. 2007; Scoville et al. 2007a, b). These 
correlations largely arise due to choices in how the image 
resampling and combination process is done. Figure 2 in 
Rowe et al. (2015) illustrates the noise correlation function 
for these noise fields; as shown, when moving one pixel in 
the horizontal or vertical direction, this correlation function 
is ~ 0.4 to 0.5, while moving two pixels away results in it 
decreasing by another factor of 10 . 

This specific noise correlation function is adopted as a 
practical example, but we caution that the details of our re¬ 
sults (for example, detection significances without vs. with 
correlated noise) depend on the adopted level of correlations, 
which will in general differ for different datasets and anal¬ 
ysis methods. We explore variations in our default level of 
noise correlation in Section 4.5 and find that our overall con¬ 
clusions are valid for significantly stronger and weaker noise 
correlations, as well. Thus, while coadded ground-based data 
will typically have weaker correlation patterns than our de¬ 
fault (since the images are typically resampled to a pixel 
scale comparable to the native one, rather than a signifi¬ 
cantly smaller one), our overall conclusions should apply to 
typical ground-based data. 

Finally, a subset of our simulations include anisotropies 
in the correlated noise. Correlated noise is in general 
anisotropic, due (for example) to direction-dependent cor¬ 
rections for charge transfer inefficiency and camera distor¬ 
tion. Our tests reveal the degree to which these low-level 
anisotropies in the correlated noise field can contaminate 
per-galaxy shear estimates. 


3 SIMULATION AND ANALYSIS METHODS 

Image generation (Sec. 3.1) and the measurement process 
(Sec. 3.2) is performed in a single seamless pipeline using 
GalSim. Galaxy and PSF profiles are convolved and ren¬ 
dered using Fourier-based rendering methods (Rowe et al. 
2015). Many independent noise realizations (up to 10^) are 
generated for each value of SNR, and analyzed using the 
shape measurement methods in Sec. 2 . 1 . 2 . 

3.1 Simulation generation 

3.1.1 Galaxy Profile 

For this work, galaxy light profiles are represented using the 
Sersic profile (Sersic 1963), a family of light profiles that can 
describe a wide range of galaxies. The surface brightness of 
a Sersic profile varies as 

(9) 

1/2 

where F is the total flux, ri /2 is the half-light radius, and 
a(n) and h{n) are known functions with numerical solutions. 
Flux normalization determines a(n), while h{n) is set by 
requiring that the half-light radius enclose half of the flux. 


n 

1,2, 3, 4, 5 

ei 

-0.45, 0.45 

62 

-0.45, 0.45 

G/2 

0.4", 0.7", 1" 

PSF 

Kolmogorov (FWHM: 0.7") 

Pixel Scale 

0.2"/ pixel 


Table 1. Parameters used for our main set of simulated images 
from ground-based telescopes. 



Figure 1. Examples of postage stamp images created by GalSim, 
for a single choice of Sersic profile parameters. The color-scale is 
linear in the flux. From left to right, the panels show the image 
for z^ideal = cxD (no noise added) and that same image with un¬ 
correlated noise, correlated noise, and correlated noise that was 
sheared with (ei = 0 . 5,62 = —0.5). All but the noise-free image 
have noise variance chosen such that i^ideal = 15 (one of the lowest 
SNR values used for this study). The case of sheared correlated 
noise is only used in Sec. 4.5.2. 

Our tests use sixty possible Sersic galaxy profiles with 
one of five values of Sersic index, two values for each shape 
component, and three values of size. Each galaxy consists of 
a single Sersic component. Table 1 lists the values used for 
these parameters. 

3.1.2 PSF and pixel response 

The simulated PSF is a circular Kolmogorov profile with a 
full width at half maximum (FHWM) of 0.7". This choice 
is representative of long-exposure times from ground-based 
telescopes at sites with good imaging conditions. Our sim¬ 
ulated ground-based images have a pixel scale of 0 . 2 ^'' per 
pixel. 

3.1.3 Image rendering 

The galaxy and PSF (including the Kolmogorov profile and 
pixel response) are convolved using Fourier-space rendering 
in GalSim, and the light profiles are sampled onto a 50 x 50 
pixel image using the chosen pixel scale. 

3 . 1.4 Noise field 

For each value of z^ideai, Eq. ( 8 ) was used to calculate a noise 
variance to be added to the initially noise-free postage stamp 
image. Many noise realizations with that variance are then 
generated, using uncorrelated and correlated noise as de¬ 
scribed in Sec. 2.3. Example images are shown in Fig. 1, with 
the sheared correlated noise case only used in Sec. 4.5.2. 

3.2 Measurement process 

Measurements of flux, size (cr), and shape (ei, 62 ) were taken 
on as many as 10 ^ different realizations of correlated and 
uncorrelated noise for the same galaxy image. In addition 
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7.5 10 15 20 I^iaeai 


^ideal 


Figure 2. The percentage of measurement failures as a function 
of i^ideal is shown for the case of uncorrelated and correlated noise, 
for one particular galaxy model. For i^ideal > 20 (not shown) the 
failure rate falls consistently below 1 per cent. 


Figure 3. Shape measurements as a function of the ideal detec¬ 
tion significance r'ldeaP for images with correlated and uncorre¬ 
lated noise. Results for the two shape components are shown in 
separate panels. Seijei is the relative bias in shape component 
z, plotted here as a percentage. Points are shown with errorbars, 
which are typically smaller than the points themselves. 


to quantifying biases in the above quantities as a function 
of the idealized detection significance z^ideai, the many noise 
realizations are also used to estimate the effective detection 
significances defined in Sec. 2 , namely z^e, i^size, influx- These 
are strictly below z^ideai (which assumes perfect information 
about the light profile) and lower for correlated noise than 
for uncorrelated noise. 

Flux, size, and shape measurements from a single galaxy 
image were considered failures and discarded if the code 
raised a fatal error, or if one of the measurement of flux 
or size fell outside a generalized 5cr range centered at the 
mean. This range is determined by independently construct¬ 
ing the probability distribution function (PDF) of the values 
of both log(size) and log(flux), and then estimating a\ow and 
(Thigh as the difference between the 50th and 16th percentile, 
and the 84th and 50th percentile, respectively for each PDF. 
Measured values that are less than 5criow below the mean, or 
more than 5(Thigh above the mean, are excluded as failures. 

Any dataset consisting of the measurements of hux, 
size, and shape for ah noise realizations of a particular 
galaxy model and noise variance with a failure rate above 
20 per cent is discarded, out of a concern that the remain¬ 
ing 80 per cent of successful measurements might be a non¬ 
representative sample. Failure rates this high only occur for 
the case of correlated noise with z^ideai = 7.5 and 10; the fail¬ 
ures arise primarily due to a failure of the adaptive moment 
routine to converge in a reasonable number of iterations at 
low effective SNR. As a result, there are two fewer points for 
correlated noise results than for uncorrelated noise on every 
plot in Section 4. Failure rates as a function of z^ideai are 
shown for one particular galaxy model in Fig. 2, which shows 
a rapid convergence to a very low failure rate at z^ideai = 20 
for both correlated and uncorrelated noise. In principle, a de¬ 
pendence of the failure rate on galaxy position angle and/or 
shape could result in a form of shear selection bias that de¬ 
pends on the type of noise (correlated vs. uncorrelated). We 
will explore this possibility in Sec. 4.4.2. 


4 RESULTS 

Our results are organized as follows. Sec. 4.1 shows the biases 
in galaxy shapes for a single galaxy model as a function 
of the ideal detection significance z^ideai for correlated and 
uncorrelated noise. Sec. 4.2 shows our results for the effective 
detection significance defined in various ways, again for both 
correlated and uncorrelated noise. 

An obvious question to ask is that since uncorrelated 
and correlated noise with different variance (zero-lag value 
of the noise correlation function) results in different effec¬ 
tive detection significance but the same z^ideai, are the shear 
biases determined in Sec. 4.1 better quantified in terms of 
the effective detection significance? And, if defined that way, 
do the results with uncorrelated and correlated noise follow 
identical trends? We address this point in Sec. 4.3. The re¬ 
sults up to that point are for a specific galaxy model and 
to a noise correlation structure similar to that seen in COS¬ 
MOS weak leasing images. The dependence of these results 
on galaxy properties such as size, Sersic index, and shape 
(Sec. 4.4), and on PSF and noise field properties (Sec. 4.5) 
are shown next. Finally, the extension of our results to an¬ 
other galaxy shape measurement code, imSshape, is shown 
in Sec. 4.6. 


4.1 Shape bias 

Unless otherwise specified- i.e., in Sec. 4.4 and 4.5 - the plots 
and data presented in this and subsequent subsections are 
for a galaxy with ( 01 , 62 ) = (0.45,-0.45), ri /2 = 0.7'^ and 
n — 2. These midrange values were chosen to represent the 
overall conclusions drawn from analysis of ah combinations 
of galaxy parameters in Table 1. 

Fig. 3 shows the fractional bias in galaxy shape mea¬ 
surements^ as a function of the ideal detection significance, 

^ For a weak lensing measurement, what actually matters is the 
bias in the ensemble average shear rather than the per-galaxy 
shapes. The bias in that ensemble shear estimate is a weighted 
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i^ideai, for this galaxy model. The results are quite similar 
for both components, so future plots of shear biases are re¬ 
stricted to a single component. There is clear convergence 
to a bias of zero at z^ideai ~ 100 for both the correlated 
and uncorrelated cases. At lower detection significance, the 
bias becomes negative, indicating an overly round shape. 
The sign of this effect is consistent with simulation-based 
investigations of the shear biases of the re-Gaussianization 
method (Mandelbaum et al. 2012, 2015). 

At a fixed value of z^ideai, the bias is stronger for corre¬ 
lated noise than for uncorrelated noise, but the curves have 
approximately the same shape. This trend suggests an inter¬ 
pretation in terms of some effective, empirically-determined 
detection significance, an idea that will be explored further 
in Sec. 4.3. 

4.2 Effective detection significance 

This section contains empirical estimates of the effective de¬ 
tection significance z^eff (influx, J^size, and z^e), using the same 
galaxy model as in the previous subsection. The results 
shown here use 10 ^ noise realizations in order to robustly 
measure uncertainties in measured quantities. 

Measurements for a range of noise variance values span¬ 
ning over two orders of magnitude (7.5 ^ z/ideai ^ 1000), 
shown in Fig. 4, reveal that for the majority of that range, 
the ratio of effective to ideal detection significance z/gff/z^ideai 
is consistent with a single constant value that is significantly 
below 1, as in Mandelbaum et al. (2015). That ratio takes 
on different values in the case of uncorrelated noise and cor¬ 
related noise, with the former being higher (as expected) by 
roughly a factor of 2. For z^ideai ^ 20, this simple z/gff oc z^ideai 
relation is violated, with z^eff falling below the linear relation. 

For fixed z/ideai > 20, the z/gff values satisfy z/flux < < 

z/gize- This hierarchy is not necessarily consistent for galaxy 
models with a different size or shape, or for different corre¬ 
lated noise fields, as will be shown in subsequent sections. 

4.3 Remapping shape biases with z^eff 

The results from Sec. 4.2 can be used to reinterpret those 
from Sec. 4.1 in terms of effective detection significances z/gff. 
Starting with the results in Fig. 4, we can scale the values 
of z/ideai for each dataset by (z^eff/^^ideai) to plot galaxy shape 
biases as a function of z^eff. This ratio, (r'eff/i^ideai), is taken to 
be a single number that can be read off from Fig. 4 for each 
type of detection significance and noise^. The key step is to 
do this rescaling separately for uncorrelated and correlated 
noise. Our goal is to determine whether, after doing so, the 

average bias of the per-galaxy shape biases, which will in general 
depend on SNR, resolution, morphology, and other quantities. 
Nonetheless, for this study we focus on galaxy shape biases, as a 
simple way to gain some initial insight into the impact of corre¬ 
lated noise. 

^ One might ask why we do not simply use the measured values 
of z^eff directly, instead of applying the simple linear scaling z/gff oc 
^ideal- We adopt this simpler procedure because in a study aimed 
at calibrating shear biases, it seems more likely that one would 
determine roughly how the detection significance scales with the 
square root of the noise variance, and apply an overall correction 
factor, just as we are doing here. 


highly offset curves in Fig. 3 for correlated vs. uncorrelated 
noise lie on top of each other. We focus on the region where 
the shape bias is significantly nonzero, z^ideai ^ 40, since 
that is the regime where noise bias must be calibrated using 
simulations. 

Fig. 5 shows the relative bias in one shape component as 
a function of the three z^eff values and z^ideai, for both uncorre¬ 
lated and correlated noise. In Fig. 3, differences in the shear 
calibration bias for uncorrelated noise vs. correlated noise 
at fixed z^ideai were typically factors of ~ 6 . It is clear that 
the use of any of our z/gff options in Fig. 5 significantly re¬ 
duces this difference to at most 20 per cent. However, of the 
three options, z/flux most reduces the difference between the 
shear biases with correlated and uncorrelated noise. When 
plotting the bias in galaxy shape as a function of z/flux, the 
curves for uncorrelated and correlated noise lie nearly on 
top of each other, z^size is the next best option, and Ue is, 
interestingly, the worst option of the three. 

The results of this subsection suggest that when ap¬ 
plying simulation-based noise bias corrections to galaxy 
shape estimates, results from noise correlation-free simu¬ 
lations could be used to correct the results in data that 
includes noise correlations, as long as the impact of noise 
correlations on the detection significance are accounted for. 
To quantify what error might be induced by doing so, we 
take the points in Fig. 5 for uncorrelated noise, and find a 
three-parameter fitting function that describes the bias as a 
function of z/eff,ucn, using z/flux to define the effective detec¬ 
tion significance. This function is 

IT Ol 1 

Bias (%) = -0.074 - ^29-. 

^effjUcn 

In reality the expansion is more likely a constant plus 
l/^eff,ucn pfos higher order terms as shown in Hirata et al. 
(2004), but we used a single term with varying power-law 
index to effectively include some higher order behavior with¬ 
out introducing many degrees of freedom into the fit. We 
evaluate this function at the values of z/eff,cn in Fig. 5, and 
calculate the residual bias for the correlated noise case, 
defining the residual as the actual bias on Fig. 5 minus the 
bias from the fit. The RMS value of the residual bias over 
all points on the curve is 0.3%, more than an order of mag¬ 
nitude below the typical biases for z^eff < 20 . 

4.4 Dependence on galaxy properties 

This section illustrates the dependence of our results on the 
properties of the simulated galaxies, going beyond the sin¬ 
gle galaxy model used in Sec. 4.1. We test how the ratio 
(i^eff,ucn)/(i^eff,cn) inferred from Fig. 4, and the efficacy of 
the remapping described in Sec. 4.3, varies with different 
properties described below. 

4 - 4 -^ Sersic index and size 

The same process as in previous sections was carried out 
for fifteen different combinations of parameters (five values 
of Sersic index, three values of ri /2 and a single value of 
( 61 , 62 ) = (0.45, —0.45)) originally listed in Table 1. 

Carrying out the exercise in Fig. 4 to find the rela¬ 
tionship between z^eff/z^ideai for various choices of z/gff reveals 
variations in those ratios that are as large as 40 per cent 
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Figure 4. Three different measures of the effective detection significance i^eff (top l^ft: i/g; lower left: r'fluxJ lower right: i^'size) are shown 
as ^'eff/^ideal as a function of i^ideal- The average values of z^eff/^ideal across all z^ideal values is shown as (z^eff/^ideal) (horizontal lines, 
with the average value given directly below each line). Only i/ei is plotted as Ue, since it is identical to results for the other shape 
component. A “Scaled r'eff,ucn” Hae was constructed by multiplying the red points (correlated noise) by the ratio of the (z^eff/^ideal) 
lines for uncorrelated vs. correlated noise. This allows for an easier comparison between the shapes of the curves for correlated and 
uncorrelated noise. 
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Figure 5. Biases in shape measurements as a function of i/eff? for 
different choices of z/eff (^'e, r'flux? and resize)^ for both uncorrelated 
and correlated noise. For comparison, the results from Fig. 3 as 
a function of z^ideal are also shown. Thus, the plot shows just two 
curves (correlated and uncorrelated noise biases), each repeated 
at four distinct sets of positions along the x axis, with horizontal 
rescalings determined from Fig. 4. Results were determined using 
10^ noise realizations. Note the more limited horizontal axis range 
compared to Fig. 3, chosen to accentuate the region of the curves 
with a nonzero bias. 


across the entire range in Sersic index or ri/ 2 - For a fixed 
value of z^ideai (i-e., noise variance) and galaxy flux, z/gff tends 
to be largest for small galaxies (low ri/ 2 ) and low Sersic in¬ 


dex. This makes sense, since the light prohle covers fewer 
pixels and thus has effectively less noise in that case. 

Tests of the impact of correlated vs. uncorrelated noise 
use the ratio (r'eff,ucn)/(i^eff,cn), as shown in Fig. 6. The left 
column shows that at hxed Sersic index, this ratio is a weakly 
increasing function of ri/ 2 , for ah three types of The 
strength of the trend is larger at low hxed n, but still less 
than 10 per cent in the most extreme case (and therefore 
weaker than the individual trends in r'eff,ucn or r'eff,cn with 
ri/ 2 ). The right column shows that at hxed size, this ratio 
is a weakly decreasing function of Sersic n, again reaching a 
maximum variation of 10 per cent across the entire range of 
galaxy parameters and being strongest for the largest galax¬ 
ies considered. 

This stability in (r'eff,ucn)/(i^eff,cn) is fairly remarkable 
given that the range of galaxy sizes is a factor of 2.5 and 
the range of Sersic n goes from 1 to 5. The stability derives 
from the fact that z^eff,ucn/i^ideai and z^eff,cn/z^ideai both have 
stronger trends with galaxy properties as mentioned above, 
but the trends are sufficiently similar that they largely can¬ 
cel out in the ratio r'eff,ucn/i^eff,cn- Any scheme that aims to 
correct for calibration biases in shape measurements using a 
combination of simulations with uncorrelated and correlated 
noise would need to account properly for the differences be¬ 
tween Z^effjUcn/^ideal and Z^effjCn/^ideal • 

The parameters of the lines shown in Fig. 6 can be found 
in Table 2. We also conhrm that the success in remapping 
shape bias vs. z^eff to match up the uncorrelated and corre¬ 
lated noise results (Sec. 4.3) carries over to ah other values of 
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n and ri /2 tested in this section. A single example is shown 
in Fig. 7. 


4.4-2 Shape Dependence 

The next test was for the dependence of these results on 
the galaxy shape using galaxies with Sersic n = 2 and 
^ 1/2 = 0.7", for ei = 0.1,0.3,0.45 and 62 = 0. Since our 
results in prior sections were very similar for the two shape 
components, only ei was varied in this section. 

The ratio z^ucn/i^cn varies very little (< 1 per cent) over 
the range of shapes considered here. Moreover, Fig. 8 demon¬ 
strates that the remapping successfully reduces the differ¬ 
ences between the shape bias for correlated and uncorre¬ 
lated noise for the three values of ei considered. This result 
demonstrates that our results for ei — 0.45 can be consid¬ 
ered as a general result, independent of shape. 

As noted in Sec. 3.2, selection effects that correlate with 
shape can cause a bias in ensemble shear estimation. If this 
bias differs for correlated and uncorrelated noise, then that 
would also be a potential issue with remapping shear biases 
for the two noise types. However, for the simulated galaxy 
and PSF profiles in this work, and the range of |e| explored 
in this section, we hnd that the failure rate is totally inde¬ 
pendent of |e| for both correlated and uncorrelated noise. 
It is also independent of the galaxy position angle for fixed 
|e|. Note that in a more complex situation with anisotropic 
PSFs, it is possible that this would no longer be the case. 


4.5 Dependence on image properties 

In this section, we explore the dependence of our results 
on the imaging properties, including the nature of the noise 
correlations and the type of PSF. 


4 . 5.1 Different scale-length of noise correlations 

Our first tests are of the dependence of the results on the 
scale-length of the noise correlations. While our original 
noise correlation function corresponded to that in the COS¬ 
MOS weak leasing science images, GalSim permits users 
to “magnify” the noise correlation function. This operation 
corresponds to preserving the variance (zero-lag correlation 
function) while changing the scale length of correlations. Our 
tests involve linearly expanding the scale-length of correla¬ 
tions by factors of /exp = 0.5, 2.0, for a galaxy with Sersic 
n = 2 , shape ( 61 , 62 ) = (0.45,-0.45), and ri/ 2 = 0 . 7 ". 

The results for r'eff,ucn/i^eff,cn are shown in Fig. 9. Since 
only the correlated noise field is being modified, these should 
be interpreted as modifications in r'eff,cn with r'eff,ucn remain¬ 
ing fixed. The ratio z^eff,ucn/i^eff,cn increases linearly with 
/exp, implying that the detection significance with correlated 
noise, z^eff,cn, is smaller for large /exp. The sign of this trend 
makes sense: for a fixed point variance, larger /exp implies 
that the noise correlations extend over larger spatial scales, 
which should decrease the effective detection significance for 
extended objects. 

Fig. 10 shows that remapping the shape biases based 
on the r'eff/i^ideai successfully matches up the shape biases 
for uncorrelated and correlated noise for all three values of 


/exp. Comparison of the panels indicates that this is a non¬ 
trivial statement, since the difference between the results for 
uncorrelated and correlated noise as a function of the point 
variance can be quite large for large values of /exp. Hence 
our results seem fairly robust to the scale-length of the noise 
correlations. 

It is possible that if the scale-length of the correlations 
became comparable to the galaxy size itself, then additional 
effects would come into play. The cases considered here are 
not in this regime (which should not be relevant for too 
many galaxies that will be used for weak leasing), and we 
defer this question to future work. 


4 . 5.2 Anisotropy in noise correlations 

The correlated noise field used throughout this work has 
a mild level of anisotropy in the noise correlations (a few 
per cent). Generally, unless explicitly accounted for^ these 
anisotropies are expected to be imprinted at some level 
in the measurements of galaxy shape, particularly at low 
detection significance. This effect is explored in this sec¬ 
tion, with a galaxy with Sersic n = 2 , ri/ 2 = 0 . 7 ", and 
( 61 , 62 ) = (0.45,-0.45), using our original correlated noise 
field as well as versions that have an additional shear ap¬ 
plied to the noise correlation function. To make the trends 
visually apparent, the noise correlations are sheared by an 
exaggerated amount, 6 noise = 0.5 (see right-most panel of 
Fig. I for an example image). 

In this case, what matters are the results with shear ap¬ 
plied to the correlated noise compared to the results with the 
original correlated noise. More specifically, when measuring 
6 z,orig for shape component i and the original noise field, 
and 6 i,sh with the sheared noise field, for the case where the 
expected shape is 6 i,exp, the following differences are calcu¬ 
lated and plotted: 

Gj^exp CijOrig Cj^exp 
Gz,exp C-j^exp 

Gj^orig 
G-ijexp 

The results are shown in Fig. II. When shearing the 
correlated noise in one component, the dominant effects on 
the shape bias are in that component, but they are not com¬ 
pletely zero for the other component. The sign of the effect 
is as expected, with positive noise shear in one component 
translating into a positive shape bias. For the lowest r'ldeai 
value considered (z^ideai = 15), the magnitude of the addi¬ 
tional shape bias is typically of order 5 per cent. We expect 
this to scale linearly with the noise shear, so for more real¬ 
istic values like 6 noise = 0 . 1 , the maximum shape bias would 
be of order I per cent for the lowest SNR galaxies in our 
simulated sample. This effect is therefore large enough to 
matter for upcoming surveys that require understanding of 
shear calibration biases at the sub-percent level. 

An additional source of bias not considered here is se¬ 
lection bias. If there is anisotropic correlated noise, and if 


^ Most shear estimation routines do not have a means to account 
for correlated noise fields, regardless of their anisotropy level. 


MNRAS 000, 000-000 (0000) 





Correlated noise 9 









ri/2 


- 

A Z/fiux 

- 

- 

▼ 

- 

- 

☆ Z^e 

- 

A - A t * * 



* * » - ▲ 




r^fo=OAO" 




Figure 6. The ratio of detection significance for uncorrelated vs. correlated noise, (^'eff,ucn}/{^eff,cn} (labeled more simply as v'ucn/i^cn), 
as a function of galaxy properties. The left (right) column show trends with galaxy size, ri/ 2 i at fixed Sersic n (with Sersic n at fixed 
^ 1 / 2 )- The equations for the plotted trendlines can be found in Table 2. 


n = 1 

d 

II 

Flux Size Shape 

Flux Size Shape 

y = 0.16r'^2 + 2-16 y = 0.25r'^2 + 2-03 y = 0.29r'^2 + 1-95 

y = -O.Oln' + 2.09 y = -O.Oln' + 1.92 y = -0.02n' + 1.83 

n = 3 

N 

d 

II 

Flux Size Shape 

Flux Size Shape 

2/= 0.09r;^2+2-11 y = 0.16r;^2 + 1-98 y = 0.17r[^^ + 1.S8 

y =-0.02n' + 2.12 y =-0.02n'+ 1.99 y =-0.03n'+ 1.89 

n = 5 

q 

II 

Flux Size Shape 

Flux Size Shape 

y = 0.06r;^2+2-09 y = 0.12r;^2 + y = 

j/ =-0.03n' + 2.15 y =-0.03n'+ 2.03 y =-0.04n'+ 1.94 


Table 2. The equations of the trendlines shown in Fig. 6, organized according to their positions in that plot. Here, y represents i^ucn/^^cn 
while n' and r '^2 represent the varying quantities respectively, offset such that n' = (n — 2) and = (^ 1/2 “ 0.7"). 


galaxies are more or less likely to be selected for measure¬ 
ment if their shape aligns in some particular way with re¬ 
spect to that noise field, then this will cause an additional 
bias in ensemble shear estimates that is not captured by our 
tests. 


4 . 5.3 Ground-like versus spaee-like data 

To test the impact of the PSF type and sampling, simu¬ 
lated space-based data was also generated, with parame¬ 
ters in Table 3, and used for the same analysis as was al¬ 
ready presented for the ground-like simulations. For a galaxy 
model with ( 61 , 62 ) = (0.45,-0.45), ri /2 = 0 . 2 ", and Sersic 
n = 2 , the ratio of r'eff/^^ideai is qualitatively similar to that in 
Fig. 4, in the sense that there are flat curves for z^ideai ^ 20. 
The actual ratios differ, with values of r'e/i^ideai = 0.83 
and 0.42 for uncorrelated and correlated noise, respectively; 


influx/z^ideai = 0.43 and 0.20; and z^size/z^ideai = 0.47 and 0.23. 
Thus, in all three cases, the ratio of z^eff for correlated vs. for 
uncorrelated noise is ~ 2. However, the hierarchy from Fig. 4 
(z^flux < z/e < z/gize) IS not rcspccted here, implying that the 
hierarchy depends on the details of the type of data. Note 
that the galaxy model here is smaller than for ground-like 
data, with ri /2 = 0.2" instead of 0.7". This choice, which is 
meant to reflect the fact that the typical galaxy that is an¬ 
alyzed in space-based datasets is smaller than can be easily 
resolved in ground-based datasets, may in part be responsi¬ 
ble for the change. 

Despite this difference in the error hierarchy. Fig. 12 
clearly shows that our conclusion about remapping shape 
biases using z/gff to account for the differences between un¬ 
correlated and correlated noise is also valid for space-like 
data. In this case, the best remapping is provided by z/gize, 
while the worst comes from i/q. 


MNRAS 000 , 000-000 (0000) 







































10 


Gurvich & Mandelbaum 



Figure 7. Biases in shape measurements as a function of i/, for 
different choices of z/eff (^e, ^flux? and z^size) and i^ideal: for both 
uncorrelated and correlated noise. Unlike Fig. 5, this plot uses a 
different galaxy model with n = 5 and ri /2 = • The remapping 

by i^eff is still successful in reducing the horizontal offset in the 
galaxy shape bias curves for correlated and uncorrelated noise. 


n 

1,2, 3, 4, 5 

ei 

-0.45, 0.45 

62 

-0.45, 0.45 

rii2 

0.1", 0.2", 0.3" 

PSF 

Airy {X/d: .041253") 

Pixel Scale 

0.03'7 pixel 


Table 3. Parameters used for simulated space-based images. 

4.6 Generalizing to other methods 

An obvious question about the results earlier in this sec¬ 
tion is how specific are they to re-Gaussianization, versus 
applying more generally. For one galaxy model, the same 
set of tests were carried out using imSshape, a maximum- 
likelihood forward model-fitting approach. The configura¬ 
tion settings used to run imSshape for this test are given 
in Appendix A. The galaxy model used for this test has the 
following parameters: e = (0.45,0.45), n = 2, ri /2 = 0.7", 
using our standard PSF and pixel scale for ground-based 
simulations from Table 1. 

First, the shape bias as a function of z^ideai is shown in 
Fig. 13. The results with uncorrelated noise can be com¬ 
pared against those from the upper left panel of figure 2 in 
Kacprzak et al. (2012), which investigated noise bias in the 
context of imSshape. While the shape of our curve appears 
different from the one in Kacprzak et al. (2012), it is impor¬ 
tant to keep in mind that our results go to lower signal-to- 
noise ratio. For the range that the plots have in common, the 
shape of the curve and magnitude of the bias at low detec¬ 
tion significance is actually fairly consistent (with some mild 
differences that could arise due to the galaxy models not be¬ 
ing identical). The curve with correlated noise does not seem 
to go to positive values, unlike the curve with uncorrelated 
noise; however, the curves have in common the behavior at 
the lowest detection significance, where they both go to neg¬ 
ative values. 

Next, we investigated the z^eff for shape, size, and flux 
from imSshape (Fig. 14). It is important to bear in mind 
that z^size for imSshape is intrinsically different from z/gize 



Figure 8. Biases in shape measurements as a function of z/, for 
different choices of z^eff (z^e, cind i^size) cind i^ideal: for both 

uncorrelated and correlated noise. Unlike Fig. 5, this plot shows 
results for two other values of galaxy shape, for which the remap¬ 
ping by z/gff is successful in reducing the offset in the galaxy shape 
bias curves for correlated and uncorrelated noise. 


from adaptive moments, since the former is intrinsic size and 
the latter is PSF-convolved. In general, z^gize from imSshape 
is therefore expected to be lower than that from adap¬ 
tive moments. Indeed, Fig. 14 reveals quite low values of 
i^size/i^ideai. It is Icss clcar why z/fiux/i^ideai IS SO low, imply¬ 
ing that detections that are visually quite apparent in images 
have a flux-based detection significance around 1. Also, it is 
generally the case that z/gff/z^ideai (for all three types of z/gff) 
are not consistent with being flat, unlike the results in Fig. 4 
from adaptive moments, though Usize is closer to flat than 
Ue- Nonetheless, our final step is to check the remapping of 
the galaxy shape biases from Fig. 13 using z/gff in Fig. 15. 
Clearly the remapping brings the curves significantly closer 
together, though their behavior is not identical due to the 
lack of a small positive bias for some range of detection sig¬ 
nificance when using correlated noise. Thus, while the results 
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Figure 9. The ratio ^'eff,ucn/^eff,cn is shown as a function of /exp, 
the linear expansion factor that was applied to the scale length 
of the correlated noise field. The lines are the result of a linear fit 
to the points. 

for remapping from re-Gaussianization do not completely 
carry over to imSshape, the remapping is still somewhat 
beneficial in reducing differences in shape biases for uncor¬ 
related and correlated noise. 


5 CONCLUSIONS 

In this work, we have investigated the impact of corre¬ 
lated noise fields on galaxy shape estimation using the 
re-Gaussianization galaxy shape estimation method. Most 
galaxy shape estimation methods ignore correlated noise, 
which may be present due to galaxies below the detection 
threshold, or image resampling and other elements of the 
image processing. Ignoring correlated noise leads to overes- 
timation of the detection significance, and could in principle 
lead to incorrect conclusions about biases in galaxy shape 
measurement. 

Our results suggest that correlated noise tends to 
change the effective detection significance in ways that de¬ 
pend only mildly on galaxy properties, and are easily charac¬ 
terized. Moreover, while the galaxy shape biases as a func¬ 
tion of the variance of the noise field can differ dramati¬ 
cally for uncorrelated and correlated noise, using an empir¬ 
ical detection significance estimator that takes into account 
the impact of noise correlations can account for essentially 
all of this difference. This statement is robust to significant 
changes in the parameters of the galaxies, PSFs, and corre¬ 
lated noise fields in the images. While this result is perhaps 
unsurprising given the underlying cause of noise bias (as 
discussed in Sec. 1), a direct validation using simulations 
is nonetheless valuable in confirming our understanding of 
noise bias. 

From a theoretical perspective, this result suggests that 
the correlated noise is not doing anything untoward to the 
images, and the primary reason they change the results is 
due to the degradation in detection significance for extended 
objects. As a practical matter, this is useful; depending on 
the required level of accuracy, it may not be necessary to run 
suites of simulations to calibrate noise bias for many differ¬ 
ent correlated noise fields, but rather that simulations with 



Figure 10. Biases in shape measurements as a function of 
for different choices of r'eff (^e, ^flux, cind i^size) cind z^ideaU for 
both uncorrelated and correlated noise. Unlike Fig. 5 (which has 
/exp = 1, i-e., no modification in the correlation length of the 
noise correlations), this plot shows results for two different values 
of /exp, a linear expansion of the scale length of noise correlations, 
for which the remapping by is successful in reducing the offset 
in the galaxy shape bias curves for correlated and uncorrelated 
noise. 

uncorrelated noise (or a single level of correlated noise) can 
be used to calibrate noise bias even when the level of noise 
correlations is not completely uniform across the dataset. 
This result rests on using a reliable estimator of detection 
significance that includes the effect of noise correlations. 
The development of such a calibration scheme that prop¬ 
erly takes all factors into account is beyond the scope of this 
work. In addition, the impact of anisotropy in the correlated 
noise may cause shear biases that are significant enough to 
matter for large upcoming weak leasing surveys when using 
galaxies near the detection limit. 

Defining the detection significance can be non-trivial, 
and if it is done in a way that correlates with shape, this 
can result in a selection bias with respect to the leasing shear 
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Figure 11. Modifications to the galaxy shape bias when shearing 
the correlated noise field by the amounts given in the legend are 
shown as a function of the ideal detection significance, z^ideal- 



Figure 12. Biases in shape measurements as a function of u, 
for different choices of i/eff (^e, ^flux? cind i^size) cind ^'ideaP for 
both uncorrelated and correlated noise. Unlike Fig. 5, this plot 
shows results for space-like simulations using the parameters in 
Table 3, for which the remapping by z/eff is also successful in 
reducing the offset in the galaxy shape bias curves for correlated 
and uncorrelated noise. 
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Figure 13. Galaxy shape bias for a single galaxy model as a 
function of z^ideal for im3shape, for uncorrelated and correlated 
noise models. 


Our results about the impact of correlated noise may 
explain some of the findings of Hoekstra et al. (2015), who 
note that the estimated shear calibration bias for galaxies 
limited to r-band magnitude < 25 depends on the limiting 
magnitude of the simulated galaxy sample, only converg¬ 
ing at a limiting magnitude fainter than 26.5. The galaxies 
below the detection limit, while not directly used for shear 
estimation, affect the results for the brighter galaxies that 
are used, presumably due to the fact that they constitute 
a source of correlated noise. If this is the primary reason 
for the need to include galaxies below the detection limit 
in the simulations in order to get the shear biases to con¬ 
verge, it is possible that the simulated galaxy sample could 
actually be truncated slightly below (rather than well be¬ 
low) the detection limit if the noise held is generated with 
a level of noise correlations consistent with the inclusion of 
fainter galaxies. Including the deblending issues that arise 
from objects right below the detection threshold will require 
the simulated sample to go at least ~ 0.5 magnitudes fainter, 
but it is worth investigating whether the rest of the effect 
could be incorporated using correlated noise. This simpli- 
hcation would make the simulation process less expensive, 
since the steepness of the galaxy number counts means that 
the sample of galaxies that must be simulated is enlarged 
by a substantial factor when going 1.5 magnitudes fainter. 

The generality of our result that shape biases can be 
remapped from uncorrelated to correlated noise is unclear, 
though our preliminary results suggest that it is benehcial 
for one other method of galaxy shape estimation and math¬ 
ematically there are reasons to expect it to be more widely 
valid. More generally, our results highlight the complexity 
of noise bias in the case of correlated noise, and the need 
to account for its impact on effective detection signihcances 
when calibrating noise bias in weak leasing measurements. 


(e.g. Jarvis et al. 2015). While we have verihed that shape 
measurement failure rates do not tend to depend on the 
magnitude of the galaxy ehipticity for either type of noise, 
it is unclear whether the correlated noise could interact in 
some way with the dehnition of detection signihcance in a 
way that modihes selection biases. That issue is beyond the 
scope of this paper. 
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Figure 14. Using im3shape measurements, three different measures of the effective detection significance i/eff (^e, i^size: and reflux) are 
shown as r^eff/^ideal as a function of z^ideal- The average values of z^eff/^ideal across all i^ideal values is shown as (z^eff/^ideal) (horizontal 
lines, with the average value given directly below each line). A “Scaled r'eff,ucn” liae was constructed by multiplying the red points 
(correlated noise) by the ratio of the (z^eff/^ideal) lines for uncorrelated vs. correlated noise. This allows for an easier comparison between 
the shapes of the curves for correlated and uncorrelated noise. 
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Parameter Value 


model_name 

sersics 

perform_pixel_integration 

N 

verbosity 

-1 

minimizer_verbosity 

-1 

noise_sigma 

1.0 

background_subtract 

N 

upsampling 

7 

n_central_pixel_upsampling 

0 

n_central_pixels_to_upsample 

0 

padding 

0 

stamp_size 

50 

rescale_stamp 

Y 

minimizer_max_iterations 

150 

levmar_epsl 

le-40 

levmar_eps2 

le-40 

levmar_eps3 

le-40 

levmar_eps4 

-1 

levmar_eps5 

le-05 

levmar_tau 

le-5 

levmar_LM_INIT_MU 

le-20 

minimizer_loops 

1 

sersics_bulge_A_min 

-50 

sersics_radius_start 

10 

psf input 

psf image single 


Table Al. Parameters used for im3shape. 


APPENDIX A: IM3SHAPE SETTINGS 

Table Al shows the parameters used for all imSshape plots 
used in this work. 
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